Morphology of ledge patterns during step flow growth of metal surfaces vicinal to 
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The morphological development of step edge patterns in the presence of meandering instability 
during step flow growth is studied by simulations and numerical integration of a continuum model. 
It is demonstrated that the kink Ehrlich-Schwoebel barrier responsible for the instability leads to an 
invariant shape of the step profiles. The step morphologies change with increasing coverage from a 
somewhat triangular shape to a more flat, invariant steady state form. The average pattern shape 
extracted from the simulations is shown to be in good agreement with that obtained from numerical 
integration of the continuum theory. 
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Epitaxial growth on vicinal surfaces is known to give 
rise to interesting growth instabilities under suitable con- 
ditions, e.g. to step bunching, mound formation and 
meandering of the step edgesB. The meandering insta- 
bility emerges when the interlayer mass transport from 
the upper side of the step is reduced due to the Ehrlich- 
Schwoebel barriern enhancing growth of protrusions at 
the step edges. This is now known as the Bales-Zangwill 
instability (BZI)El which tends to destabilize the ledge 
morphology due to terrace diffusion and asymmetric in- 
terlayer crossing. There is no diffusion along the ledges 
in BZI. However, recently it was found that in the case 
of 1 -|- 1 dimensional growth there is an analogous phe- 
nomenon due to the kink Ehrlich-Schwoebel harrier for 
going around a kink site at the step edge. The corre- 
sponding kink Ehrlich-Schwoebel effect (KESE) leads to 
growth of unstable structures at the step edges with a 
dynamically selected wavelengthcl. The ledge instabilities 
were originally found and-reported experimentally on the 
Cu(l,l,17) vicinal surfaccD but attributed to the BZI sce- 
nario. More recent STM experiments on the Cu(l,l,17) 
surface proposed that th* formation of the regular pat- 
terns is due to the KESEU. Since then theoretical studies 
of the meandering instability have shown that the KESE 
indeed supersed^the BZI in the formation of the pe- 
riodic patternsuij'u and eventuallyp. leads to an in-phase 
motion of the step edge structuresQ^El. 

Instability and wavelength selection of the step edge 
patterns due to the KESE have beenj-studied within the 
framework of a contitLu.um step modelcl, the solid-on-solid 
(SOS) lattice modeffl, and semi-realistic Monte Carlo 
(MC) simulationaJ. In particular, in the recent MC workQ 
it was shown that on vicinal Cu(1,1,to) surfaces the ob- 
served instability is due to the KESE and the competing 
BZI is of no importance in the length and time scales 
considered. The role of dimer nucleation in determining 
the selected wavelength was confirmed, in good agree- 



ment with the theoretical scaling relation^ and more re- 
cent SOS simulationds. In the MC simulations, there was 
evidence of phase locking of the ledge structures at the 
largest coverages studied, but this was not quantitatively 
confirmed. 

Regarding the step morphologies, a triangular shape 
has been predicted to occur for a strong KESE and a 
rounded, more flat shape far a weak KESEo. However, 
the MC simulations of Ref.Q indicate that in the case of 
a strong KESE there is in fact an interesting shape tran- 
sition from narrow, somewhat triangular shapes in the 
initial stage of growth to more rounded patterns in the 
large, coverage regime. Moreover, the MC simulations of 
Ref.El and SOS model results of Ref.El are in disagreement 
with asymptotic evolution of the step profiles as predicted 
by the continuum theoricaHllid. In this work we study 
the ledge morphologies of growing steps on the vicinal 
Cu(1,1,to) surfaces in detail. In particular, we study the 
onset of the in-phase growth and the phase-locking of 
the step profiles in the presence of a strong KESE. Our 
results show that the ledge morphologies assume an in- 
variant shape due to an interplay between various mass 
transport currents and phase-locking of the steps. We 
show how the ledge profiles from the MC simulations can 
be reproduced by explicitly including the relevant mass 
transport currents on the surface. On continuum level 
this indicates a delicate balance between the various cur- 
rents that determines and stabilizes the invariant ledge 
shapes. _ 

The model system used here is as in Ref.Q, based on 
MC simulations of a lattice gas model with energetics 
fromJtiie effective medium theory (for more details, see 
Refs.DO). Our MC method is efficient enough to sim- 
ulate growth of Cu up to ten monolayers (ML) under 
realistic temperature and flux conditions. The tempera- 
ture range explored here was T — 240 — 310 K and the 
flux F = 3 X 10^3 - 1.0 ML/s. Thus the ratio between 
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Fig. 1: Snapshots of typical ledge profiles with step ori- 
entations in the close packed [110] direction at T = 300 
K with F = 8 X 10-2 ML/s for coverages = 0.4,2.0, 
and 10.0, in figures (a)-(c) (lateral and vertical scales are 
1000 and 70 lattice spacings, respectively). In (d) the 
shape transition is shown. The profiles have been ob- 
tained by averaging over the meander periods at 6 = 0.4 
ML (circles) and at 6* = 2.0 ML (squares). The horizon- 
tal direction is scaled with the wavelength A = 120a and 
in the vertical direction with the roughness w = 1.4a and 
w = 5.6a for coverages 0.4 ML and 2.0 ML, respectively. 



a) 



0.1 





— $ — § § I 



log^^O) 

Fig. 2: (a) Lateral moments M„ of the step meander pe- 
riods (n — 2, 4, 6, 8, 10 from top to bottom) as a function 
of coverage in a semi-logarithmic scale. The moments ap- 
proach constant values already before 2.0 ML. (b) Width 
w of the profiles as a function of the coverage. No sat- 
uration is observed up to 10 ML. The slope of the solid 
line corresponds to /3 = 1/3. 



the terrace diffusion and the flux D/F ^ 6 x lO'"^ — 9 x 10'' 
in units of the lattice constant a — 0.361nm, corpsr 
sponding to a typical molecular beam epitaxy regimella. 
The energetics of the model also specifies the impor- 
tant length scales controlling step flow growth. These 
are Ic^ the length scale for dimer nucleation at the step 
edgald, and the kink Schwoebel lengthQ is = exp[(i?s — 
Ed)/kBT] — 1 which is related to the energy barriers 
Eg = 0.52 eV and Ed = 0.26 eV for jumps around a 
kink site and along a straight edge, respectively. For the 
close packed [110] ledges, £s — lO'* and £c ^ 10^ ^FfflW 
room tgmperature corresponding to strong KESEo^BO. 
In ReLQ it was shown that the wavelength of the step 
edge patterns is given by ic = {^'^E>s/ FL)" , Dg being the 
adatom diffusion constant along the straight edge, with 
a scaling exponent a w 0.23, and an effective barrier of 
i^off = 75± 10 meV. Both are in gQod agreement with the 
exact values which give a = 1/4113 and Eps — Ed/4: = 65 
meV, respectively. Our previous studylll was done on a 
Cu(l,l,17) surface but we have checked the results also 
with smaller terrace widths. 

Simulation results for the step edge profiles on 
Cu(l,l,17) are shown in Figs. l(a)-l(c) after deposition 
of 6* = 0.4, 2.0, and 10.0 ML, respectively, at T = 300 K 
with F = 6 X IQ-^ ML/s. In the beginning of growth 
(Fig. 1(a)) the shape of the patterns is somewliai trian- 
gular as predicted for a relatively strong KESEllil. The 



meandering structures are not yet completely in the same 
phase indicating that the diffusion field has not yet cou- 
pled the subsequent step edge tmins, a typical feature 
for KESE dominated meanderingS. However, a selection 
of the relatively well-defined wavelength^for all ledges is 
apparent already at this stage of growthQ. At larger cov- 
erages the meandering of steps begins gradually to phase- 
lock, seen in Fig. 1(b), and in-phase growth and phase- 
locking seem complete at largest studied coverage of 10 
ML shown in Fig 1(c). However, now the average shape 
of the patterns is clearly different from that at low cover- 
ages. The shape of the average patterns is more rounded, 
as predicted for a weak KESE. In Fig. 1(d) we show this 
change by comparing-average ledge profiles after 0.4 and 
2.0 ML, respectivelyM 

From Fig. 1 it is clear that there is no coarsening of 
the structures when the coverage is large enough. The 
steady state pattern shape seems to be governed by ge- 
ometric constraints which is a sign of asymmetry of the 
growth rates between bottom and top parts of the steps. 
This asymmetry is a general feature in .niaiiy models of 
step growth with or without coarseningtStZI. Moreover, 
a quantitative inspection of the patterns at larger cov- 
erages suggest that the profiles have an invariant shape. 
This can be seen by examining the nth lateral moments 
Mn{9) = {Ci{x,9)x^)i^x of the meander periods Q{x,6), 
shown in Fig. 2(a). The scaled even moments approach 
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Fig. 3: The average shape of the step patterns at 
T = 300 K and 9 = 8.0 ML with F = 6x 10'^ ML/s 
(squares). Changing the temperature or the flux as de- 
scribed in the text does not have any effects within the 
error bars. The sohd line is the stationary profile ob- 
tained by integration of Eq. (||). Good agreement be- 
tween the average and the integrated profiles is evident. 
The inset displays the relative differences of the profiles, 
A = (Call - Ci)Kaiu where C,aU is the profile with aU 
currents included, and i — k,e, SB denotes the solution 
with only a small contribution for the KESE current from 
the Gibbs-Thomson and symmetry breaking currents, re- 
spectively (from top to bottom in the inset). 

their steady state values already at 6* w 2 ML. The pat- 
terns are still changing, however, which can be seen from 
the roughness of the step u;(6') = iCi^: ^)'^)x^ where 
({x, 9) is the step profile. It does not show any sign of 
saturation up to the largest coverage in the simulations. 
Instead the roughness follows w{9) ~ 9^ , with [3 « 0.3 
as shown in Fig. 2(b). It is interesting to note that al- 
though the roughness does not saturate the shape of the 
periodic structures attains an invariant form. 

Our simulation results show that the profile shape is 
rather insensitive to deposition and temperature con- 
ditions. This suggests that the invariant shape is not 
dependent on the relative magnitudes of the various 
diffusion processes but rather is a result of geometric 
constraints due to crowding and in-phase evolution of 
the step edges. In order to justify this assumption 
we compare the MC profiles with continuum profiles 
which are obtained as stationary solutions to the dy- 
namic equation dtC = —dxJtot, where Jtot is the total 
mass current at the step edge. The most important 
partial currents which we take into account in the to- 
tal current here, when expressed in terms of the variable 
m(x) — {dxC)/y/^ + {dxCy and appropriately scaled, aije 
the mass current due to the destabilizing strong KESEq 
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Fig. 4: The mass currents Eqs. (l)-(4) are shown using 
the integrated profile as an input. Note the difference 
between the vertical scales for (c) and (d). 

the stabilizing curpent||d«e to the Gibbs-Thomson effect 
and edge diffusionlliHl3't3 

Je = ^^(yT^+^ydxxm)VT^, (2) 
the, nurrpnt into the step edge from the deposition 

Jd = Lm\/ 1 — m^, (3) 



and the front-back symmetry breaking current! 




J SB = ^ — (dxm)(dxxm)\/ \ - -m? 

^—m(dxm){i-Tr?). (4) 

In these expressions L is the terrace width, Ds is the 
macroscopic diffusion constant on the terrace, I3l is the 
macroscopic diffusion constant-akjng the step edge, and F 
is the step stiffness (see Refs.oEj for the definitions and 
experimental values of the parameters, respectively). All 
length scales are given in the units of the lattice constant. 
By requiring the condition of stationarity 

Jtot = Jk + Je + Jd + J SB = 0, (5) 

we obtain a second order differential equation for m{x). 
The stationary profiles are obtained by solving Eq. (|^) 
numerically for given initial conditions m(±l) = ±mo. 

The stationary solution is found using the value toq ~ 
0.97 as the boundary condition in order to match the 
end points with the slopes of the patterns obtained from 
the MC simulationscil. The other parameter values of 
the integration are based on known energetics of Cu, 
yielding = 700 - 1600, is = 2 x 10^ - 2 x 10^ 
Dl/{DsL) = 250 - 700, and DsT/F = 0.5 - 4000 in 
the range T = 240 - 300K and = 3 x 10"^ - 10"^ 
ML/s. For the step stiffness we used the expression 
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e^[EklkBT]/2, where Ek = 0.13 eV is the kink 
3. In aU cases we set L — 10 for the terrace 
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width. The resulting profiles are shown in Fig. 3 with 
various values of the parameters. The shape is rather 
independent of the details of the currents in agreement 
with simulations. In Fig. 3 the average shapes obtained 
from the simulations are plotted with a few different flux 
rates. In the inset we show how the resulting profile 
deviates from the complete one when each of the mass 
currents is forced to be small. 

In Fig. 4 we show the mass currents using the inte- 
grated profile as an input. It is now seen that for the 
invariant profile there is a delicate compensation of the 
currents, the Gibbs-Thomson current compensated by 
the sum of the KESE, the deposition, and the symmetry 
breaking currents. This compensation happens for the 
specific shape of the profile, and cannot take place e.g. 
in the case of a triangular shaped profile as obtained in 
the initial stages of growth. In determination of the sta- 
tionary profile shape the front-back symmetry breaking 
and the geometric constraints contained implicitly in the 
initial conditions are crucial. p 

In summary, the MCQ and the SOSH simulations have 
proven that the KESE is the dominant mechanism behind 
the meandering instability and that it leads to the selec- 
tion of the dominant wavelength determined by dimer 
nuclcation at step edges. In this work we have shown 
that the KESE also induces an invariant shape of the 
step profiles during in-phase growth. This occurs even 
though the overall roughness of the step structures wiff) 
shows no signs of saturation. The value of the corre- 
sponding scaling exponoat [3 ~ 0.3 is consistent with the 
case of an isolated stcpa. The SOSjnodel gives for the 
strong KESE an exponent (3 « 0.57Q, while for a-colleCj- 
tion of steps in the phase-locking regime P = 1 /SO'EJO. 
This puzzling behavior of dynamical scaling is apparently 
related to the strict in-phase growth and consequent for- 
mation of the invariant shape of the profile. The fact 
that the shape remains invariant although the roughness 
does not show any sign of saturation indicates a subtle 
coupling of the step edge currents with the stationary 
morphology. By numerically integrating the continuum 
equation we have shown how the interplay between vari- 
ous surface currents determines the invariant step shapes. 
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